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Abstract 

The quality of two different separable expansion methods (W matrix and 
Ernst-Shakin-Thaler) is investigated. We compare the triton binding energies 
and components of the triton wave functions obtained in this way with the 
results of a direct two-dimensional treatment. The Paris, Bonn A and Bonn 
B potentials are employed as underlying two-body interactions, their total 
angular momenta being incorporated up to j < 2. It is found that the most 
accurate results based on the Ernst-Shakin-Thaler method agree within 1.5% 
or better with the two-dimensional calculations, whereas the results for the 
W-matrix representation are less accurate. 

PACS number(s): 21.60.-n, 27.10.+h 
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I. INTRODUCTION 



Separable approximations or expansions of the underlying two-body interaction play 
an essential role in calculations of few-nucleon systems. In the three-body problem such an 
input reduces the two-dimensional Faddeev equations, or the corresponding Alt-Grassberger- 
Sandhas (AGS) equations, to one-dimensional effective two-body equations ffl. Analogously, 
the four-body AGS equations go over into effective three-body, and, after repeated applica- 
tion of separable expansions, into effective two-body equations || . A considerable reduction 
of the complexity of the original problem is achieved in this way. What remains to be done, 
however, is a careful test of accuracy of the respective separable representations. 

In early calculations Yamaguchi-type separable potentials were employed. Though not 
very realistic, they simulated characteristic aspects of the nuclear interaction, leading thus 
to qualitatively correct cross sections for three-nucleon scattering. Stimulated by this experi- 
ence, rather sophisticated separable approximation or expansion techniques were developed 
and have been applied successfully to realistic interactions. 

The Ernst-Shakin-Thaler (EST) expansion || is considered as a particularly powerful 
separable approach. Its efficiency has been thoroughly investigated for the Paris potential 
H). Indeed, based on high-rank separable expansions of this type |5||J, rather accurate 
predictions for the neutron-deuteron (n-d) cross section and some polarization observables 
were achieved |||7|. Later on, when direct solutions of the two-dimensional three-body 
equations became available |§, it was established that the results obtained in this way and 
with the EST approach are in excellent agreement [p|,|T0| . 



The W-matrix method |TlJ provides quite directly a rank-one separable representation of 
the two-body T matrix, which preserves the analytical properties of the original T matrix. 
For the Malfliet-Tjon (MT I+III) potential three-body bound-state and scattering results 
12 , 13] , including break-up |14[], were obtained on this basis in complete agreement with 



alternative treatments. Calculations for the Paris potential were also found to be rather 



reliable for most, but not all observables 13 
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For further applications of separable approximations or expansions, e.g., in processes as 



the photodisintegration of three-nucleon systems [15,16|, the corresponding radiative capture 



||17|| , or the pion absorption on such systems [lq|, additional tests of accuracy of the three- 
body wave functions involved appear desirable. In case of a 5 channel calculation and the 
EST method such a test has been made already in Ref. (TP| for the Paris potential, and, 
with a modified expansion method, in Ref. |2(| for the Argonne 14 potential restricted to 
j<2. 

In what follows we compare triton bound-state calculations performed by means of 
the EST expansion and the W-matrix approximation with direct solutions of the two- 
dimensional homogeneous three-body equations. In view of the sensitivity of the bound-state 
problem to differences in the two-body input, this comparison should provide a particularly 
relevant test of accuracy of these methods. 

The paper is organized as follows: Sections and |T| briefly describe the W-matrix and 
EST methods, respectively. In Sec. [TV] we give the relevant equations for the bound-state 
calculations using separable potentials (details of the direct treatment can be found in Refs. 
2"T|,|2"2"|]). Our results and conclusions are presented in the last section. 



II. W-MATRIX REPRESENTATION 

The W-matrix method [fTl| , |l3|| leads to a quite appropriate splitting of the two-body T 
matrix into one dominant separable part and a small non-separable remainder 

Tuip,p'; E + iO) = s T u ,(p,p'; E + iO) + R T u ,(p,p'; E + iO). (2.1) 

Similarly to the T matrix, the W matrix is defined by an equation of Lippmann-Schwinger 
type, however, with a modified non-singular kernel. After partial wave decomposition this 
equation is given by (we use units (he) 2 = 2/i = 1 and thus E = p 2 ) 

W^p,p';E) = U? t (p,p>) + E J dqq 2 ^ q) ~ ^ *° q 21 ' W^q,p'; E), (2.2) 
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where k is subject to the constraint k = \[E~ for E > and is kept arbitrary for E < 0. 
Here / and I are the orbital angular momenta, and 77 = (s,j;t) stands for the spin, the 
angular momentum j [ with the coupling sequence (I, s)j], and the isospin t of the two-body 
subsystem. The input entering Eq. ( |2.2j ) is related to the two-body potential matrix Vj(p, p') 
according to U^(p,p') = p~ l V^(p,p') p n . The separable part of the W-matrix representation 
([2.1D of the two-body T matrix is given by 



s T^p,p'; E + zO) = J2p l Wj(p, k; E) Al(E + iO) W^fc/, fc; £)/ 
and the nonseparable remainder reads 



(2.3) 



Rrpri 



T? l ,(p,p' ] E + tO)=p l 



WUP,P'; E)-Y: W3(p, fc; fc; E))- 1 W?(p', k; E) 



P 



The propagator A«, is given by 

t" 

where F~ n (E + iO) is a generalization of the Jost function 



(2.4) 



(2.5) 



{FUE + l V)Y l = 5 n , - dqq'q 21 



W n ,(q,k;E) 



(2.6) 



E + iO-q 2 

On-the-energy-shell and half-off-shell the separable part (|2.3| ) of Q2.1|) is identical with 
the exact T matrix. In fact, when inserting the momentum k = \fE for one or both of the 



momenta p or p' , we easily see that the remainder ( |2.4| ) vanishes. Therefore, the separable 
part of the T matrix has the same pole and cut structure as the full T matrix. This 
suggests to approximate the T matrices, entering the kernel of Faddeev-type equations, by 



the separable expression ( |2.3| ). To optimize this approximation, i.e., to minimize the effect 
of the neglected remainder, two criteria for the choice of the functional form of the free 
parameter k were developed [13|]. In the present work we apply the criterion based on the 
Schmidt norm of the kernel of the three-body equations. 

In the past this method has been used in elastic n-d scattering calculations |T2],|n 
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the breakup case [14], and more recently in the photodisintegration of the triton |T5| 



III. EST METHOD 



The Ernst- Shakin-Thaler (EST) method H] allows one to generate separable represen- 
tations of arbitrary rank iV that agree exactly (on- and half-off- shell) with the original T 
matrix at iV specific, appropriately chosen energies. 

For a brief outline of this method let us begin with the (partial-wave projected) 
Lippmann-Schwinger equation for the wave function, 

\^ E ) = \k E )+G (E)V\^ E ), (3.1) 

where \k E ) is the incoming wave, and Gq(E) the two-body Green's function. The dependence 
on the orbital angular momenta /, /' and on the conserved quantum numbers (the spin, the 
angular momentum j [ with the coupling sequence (I, s)j], and the isospin t of the two-body 
subsystem) is suppressed for convenience. For proper scattering solutions (on-shell), k E and 
E are related by E = k E . 

According to the EST method, a rank- A separable representation of a potential V is 
given by the form 

v est = j2 VlfaJ A,* (MV, (3.2) 

where E^, (/i = 1,...,N), is a freely choosable, but fixed set of energies. The coupling 
strengths A^ u are determined by the condition 

N 

£A^jjV|Vfc,> = V (3.3) 

Note that the "form factors" in the separable potential ( p.2| ) consist of the objects V\ip E ^), 
where \ip EfJ _) are solutions of Eq. fl3.1|) at the energies E = E^. Thus, Eq. ( |3~3"D together 
with Eq. ( |3.2|) implies that the following relation holds at the N energies E^ 

V EST \^} = V\^) = T(E,)\k E J = T EST (^)|^). (3.4) 

Here T and T EST are the two-body T matrices for the potential V and its separable rep- 
resentation V EST , respectively. Evidently Eq. (|3.4j ) means that the on-shell, as well as the 
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half-off-shell T matrices for both interactions, V and V EST , are exactly the same at the 
energies E^. 

With the form factors \g%l) = V^lip^Jji (77 = (s,j;t)) and the propagators of this 
representation, the two-body T matrix reads 

T? V {E + ^0)=Y. \ 9 y) A^(E + iO) (gll'\ (3.5) 

where 

A"(£ + iO) = ((A'')- 1 - £ (£ + 20)) _1 (3.6) 

and 

(0o(£ + zO))^ = £<<^|G (£ + iQ)\gJ)). (3.7) 
j 

For more details of this construction we refer to Refs. The chosen approximation 

energies E^ for the interaction models considered in the present study are summarized 
in Table |. These energies completely specify the separable expansion (3.2) Note that, for 
reasons of convenience, we have represented the (numerically given) form factors analytically 
(cf. Eqs. (2.1) and (2.2) of Ref. 0), and have performed the actual calculations with these 
expressions. The corresponding parametrizations can be obtained from one of the authors 
(J.H.) on request. 

In Sec. V we are going to present three-nucleon bound state results, based on the EST 
representation, with a varying rank in various two-body partial waves. We characterize these 
representations by (nin 2 n 3 ...), where n\, n 2 , n 3 , ... stand for the ranks in the 1 s , 3 Si — 3 di, 
1 Pii 3 Po, 3 Pi, 1 d2, 3 <^2, and 3 p 2 — 3 / 2 NN partial waves. For a specific rank n M in a certain 
partial wave the approximation energies E^ can be read off from Table |. They are given 
by the first entries. The notation in the results of the three-body calculations is done in 
the same order of the partial waves, but the ranks of the unused partial waves are left out. 
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IV. 3N BOUND STATE CALCULATION 



The triton bound state |^ t ) is determined by the eigenvalue equation 

(E t -H)\%) = 0, (4.1) 

3 

where the total Hamiltonian H is given by H = H + V = H + V 1 . Here we have used 

7=1 

the complementary notation V 1 = V a p for the two-body potentials, while H denotes the free 
three-body Hamiltonian. When introducing the corresponding resolvent G$(z) = (z — Ho)^ 1 
or the channel resolvents G y (z) = (z — H — Vj)" 1 , Eq. ( |4.1|) can be written in form of 
homogeneous integral equations, 

|*t> = G (E t ) V \%) = G (E t ) J2 V-y |* t > (4.2) 

7 

|* t ) = G,{E t ) V, \%) = Gj(E t ) ^(1 - 5 lP ) V p |* t ), (4.3) 



with Vy — V — V 1 being the channel interaction between particle 7 and the (af3) subsystem. 
The latter equation can also be understood as a representation of the bound state by the 
"form-factors" |F 7 ) = Vy \%), 

\%) = Gj(E t ) |F 7 ). (4.4) 

Multiplying this representation with (1 — Sp^Vy, using the relation V- 1 G 1 = T 7 Go, and 
summing over 7, we obtain for \Fp) the coupled set of homogeneous integral equations 

\F p ) = £(1 - 5 M ) T,(E t ) G (E t ) \F 7 ). (4.5) 

7 

Note that this relation may alternatively be derived by going to the bound-state poles of 
the AGS equations, providing thus their homogeneous version [I]. From (|4.2j ) and ( |4.4| ) we 
infer that the solutions of Eq. ( f4.5|) provide |\l/t) according to 

|*t> = E G o(^t) Try(E t ) G (E t ) \F y ) = E ^7) • ( 4 -6) 

7 7 

The are the standard Faddeev components, as seen by using the definition of |F 7 ) and 
the relation ( |4.3|) , 
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|V 7 ) = Go T 7 Go |F 7 ) = Go V 7 G 7 |F 7 > = G V 7 G 7 V 7 |* t ) = G K 7 |* t )- (4.7) 



Equation ( |4.5| ) will be treated numerically in momentum space, employing a complete set 
of partial- waves states \pqlbT I) . The label b denotes the set (r)KL) of quantum numbers, 
where K and L are the channel spin of the three nucleons [with the coupling sequence 
{J, \)K] and the relative angular momentum between the two-body subsystem and the third 
particle, respectively. V is the total angular momentum following from the coupling sequence 
(K, L)T, and I is the total isospin. These states satisfy the completeness relation 

oo oo 

1=E/ J dpp 2 dqq 2 \pqlbT I)(pqlbT I\. (4.8) 



bin 







The required antisymmetry under permutation of two particles in the subsystem can be 
achieved by choosing only those states which satisfy the condition (— y +s+t = — 1 . Table |I| 
contains the quantum numbers of the corresponding channels taken into account. 
Inserting the separable T matrix (|3.5| ) for the EST potentials and defining 



<(<?) = E / dpp 2 9l(p) (pqlbn\G \Fp), (4.9) 
' o 

Eq. ( |4.5|) goes over into 

oo 

F" b (q) = EE / dq'q >2 AV%(q,q',E t ) A« p (E t - \q' 2 ) F" 6 \q') , (4.10) 

with 

oo oo 

AV%(q, q', ^t) = 2 E / / dpp 2 dp' p' 2 g^p) (pqlbT I\G (E t ) \p> q' I' b' T I) g^p') (4.11) 
«' o o 

being the so-called effective potential. The recoupling coefficients entering this equation can 
be found in Ref. [|13j (or in a more compact form in [^1|] for another coupling sequence which 



can easily be changed to the present one). In case of the W-matrix representation Eqs. 
( |4.10| ) and Q4.ll ) are of similar form and, therefore, not given here. 



After discretization Eq. ( |4.10|) can be treated as a linear eigenvalue problem, where 



the energy is considered as a parameter which is varied until the corresponding eigenvalue 
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equals unity. The eigenvalues can be found by using standard numerical algorithms. A better 
approach is an iterative treatment, known as "power method" |[23|| , which is justified due 
to the compactness of the kernel of the integral equation employed. It was found that this 
method is much faster than standard eigenvalue algorithms and yields the same accuracy. 
For the direct solution of the two-dimensional Faddeev equations we used a Lanczos-type 



algorithm |||,[25[] that is even more efficient. 



For the integration in Eq. fl4.10|) a standard Gauss-Legendre mesh was chosen. The 



angular integration in the effective potential was done with 16 grid points. Table |IJ contains 
results for the Paris (EST) potential for different numbers of partial wave and an increasing 
number of mesh points for the q integration. In all cases 36 grid points were sufficient to 
get the binding energy up to 5 significant figures. For all further calculations we have used 
40 mesh points to get wave functions of high accuracy. In the calculations we also included 
q = to avoid extrapolations for small momenta in further applications. The same was done 
for both variables in the calculation of the wave function. The binding energies obtained 
in EST and W-matrix approximation for different potentials are given in Tables [TV]-|VI|, 
compared with results from a two-dimensional treatment of the Faddeev equations. 



The whole wave function can now be calculated by either using Eq. (f4.6p, or by applying 



the permutation operator P on one Faddeev component |2T 



|* t ) = (l + P)|Vi), (4.12) 

where P represents the sum of all cyclical and anticyclical permutations of the nucleons. 
From the practical point of view the latter method has to be preferred. Once is computed 
from — Go(E t ) T\{E t ) Go(E t )\Fi) the calculation of the full wave wave function via Eq. 
( [4.12| ) is independent of the rank of the separable approximation, which considerably reduces 



the computing time. 

The wave function is normalized according to 



(tt t |tt t > = £> 7 |*t> = 3 (V>i|* t > = 1- (4-13) 

7 
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It should be noted that, inserting the completeness relation (|4.8|) into (\I/ t |\I/t), one has to 
deal with an infinite number of states due to the resulting recouplings. In contrast, when 
inserting the completeness relation in (^jj^t), one has to deal with a finite number of partial 
waves corresponding to the ones in ("01 1- Plots of the wave functions for the Paris (EST) 
and Bonn A (EST) potentials are shown in Figures 1-6. The figures for the Bonn B (EST) 
potential are not distinguishable by eye from the ones for the Bonn A (EST) potential, and 
are therefore not shown. 



A. Properties of the wave function 

It is common to investigate the properties of the wave function in the LS'-coupling scheme 
(for simplicity we skip the dependence on the isospins in the notation, since they are not 
recoupled) 

\ P q{{lL)C{sS)S)TM r ). (4.14) 

In this scheme first the two orbital angular momenta and the two spins are coupled sepa- 
rately. The total orbital angular momentum is then coupled with the total spin to the total 
angular momentum of the three-body system. The total angular momentum of the triton 
is T = 1/2, the total spin S of three particles can be S = 1/2 or S = 3/2. The total orbital 
angular momentum is, therefore, restricted to C = 0, 1, 2. 

The transformation from the channel spin into the LS coupling scheme is given by 

(((lL)C,(sS)S)TM r \ 

= Y.H l+s+L+S+C+S+1 J^k T S K r S K \ ((((ls)jS)KL)TM r \, (4.15) 
ok \S j sj [T L C] 

with the abbreviation j = y / 2j~+T used only in this equation. It should be noted that for 
this transformation the wave function ( |4.12| ) has to be projected on all states that give a 
contribution due to Eq. ( [4.151 ). Otherwise the normalization constant of the wave function 



is changed. It is not sufficient to use only those channels used in the calculation of the 
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Faddeev component. Here again the recoupling of channels, as in the calculation of the 
normalization constant, plays a role. 

In the LS coupling scheme the wave function can be classified according to the contri- 
butions of the states belonging to C = 0,1,2. 

1 = 3(^1%) 

oo oo 

= 3EEE / Jdpp 2 dqq 2 {^\pq ((IL)C(sS)S)T) (pq ((lL)C(Ss)S)T\* t ) 

C S ILsS Q q 

= ]T P(£) = P(S) + P(S') + P(P) + P(D). (4.16) 
c 

The contributions to the normalization constant for a certain total angular momentum 
are denoted by P(£). In case of £ = also the symmetric and mixed symmetric spatial 
contributions P(S) and P(S'), are extracted. The antisymmetric part P(S") of the wave 
function is negligible and, therefore, has been omitted. The main contribution to the mixed 
symmetric part stems from the difference between the 1 sq and the 3 Si interaction. 



V. RESULTS AND CONCLUSIONS 

As a test of accuracy of the EST and W-matrix approaches we compare the triton binding 
energies and wave functions obtained in this way with the results of a two-dimensional 
treatment of the Faddeev equations. In Tables |V| - |V| our binding energies obtained for the 



Paris, Bonn A, and Bonn B potentials are given for different combinations of partial waves. 
For all three potentials considered here, the EST results converge with increasing rank and 
agree within 0.2% with the two-dimensional results. In case of the W matrix, the free 
parameter k was chosen differently in each partial wave according to the criterion employed 
in Ref. which consists in providing a binding energy close to the results obtained with 
other methods. This optimization method is not fully satisfactory in the j < 2 calculations, 
due to ambiguities in fixing k when including the 3 P2~ 3 f2 partial wave. The agreement of the 
W-matrix calculations with the two-dimensional results is less good than in the EST case. 
There are differences between < 4%. 
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Our five-channel, i.e., j < 1 + , calculations for the Paris (EST) potential are also in 



perfect agreement with the results by Parke et al. ||19|| . We can, moreover, compare with the 
results by Friar et al. |26| generated in coordinate space, and find again good agreement. 



For j < 2 and the Paris potential our results are exactly the same as those of Ref. [10 



The binding energies for Bonn A and j < 2 differ only slightly from those by Fonseca and 



Lehman 27 



For the Paris (EST), Bonn A (EST), Bonn B (EST) potentials the components of the 
three-body wave function, defined in the previous section, differ from the two-dimensional 
results by 0-2.5%, 0-0.05%, and 0-3% respectively. It should be emphasized, however, that 
only P(P) shows the large deviation of 3% quoted, while all other components are in much 
better agreement. 

The W-matrix results for the components of the wave function differ by 0-12%, 0-10%, 
and 0-10% in case of the Paris, Bonn A, and Bonn B potentials, respectively. Here the 
largest deviations are found both in the P(P) and in the P(S') components. 

Another test is given by the norm squared of the differences of the triton wave functions 
obtained in the separable and the two-dimensional treatments, AiV = \ \^ S ep — ^2d|| 2 - For 
the EST potentials AiV is of the order 10~ 6 , while for the W-matrix representation it is of 
the order 10~ 5 . Also in this respect the EST method, hence, leads to better results. 

Thus, we have demonstrated the high quality of the EST expansion method. For j < 1 + 



and the Paris (EST) interaction this has been done already in []19 |, and for the Argonne 
potential (with a modified expansion scheme) in [^0J. Here we have extended these former 
investigations up to j < 2, using moreover two versions of the Bonn potential. The accuracy 
achieved within the W-matrix approach is less satisfactory, namely of the order 2-10%. But 
it should be recalled that this treatment is based on a rather simple rank-one approximation 
only. 
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TABLES 



Paris potential 



partial wave 
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3 Sl - 3 d X 
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(125,2) 
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10 
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3 Po 
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(-50,2) 




l Pi 3 Po, 3 Pi 


10 


-50 


150 


300 


-150 






3 P2 ~ 3 f2 


(10,1) 


(10,3) 


(75,1) 


(75,3) 


(150,1) 


(150,3) 


(200,1) 


1 d 2 3 d 2 


10 


-50 


150 


300 


-150 







TABLE I. Approximation energies E^ used in the EST representations of the Paris, Bonn A 
and Bonn B potentials, refers to the deuteron binding energy. are lab energies in MeV. In 
case of coupled partial waves, the boundary condition chosen for the angular momentum L of the 
initial state (cf. Ref. Q) is also specified. 
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Channel Subsystem I s j n r K L 



1 


'so 








0+ 


1 


1/2 





2 


3 *i 





1 


1+ 





1/2 





3 


3 *i 





1 


1+ 





3/2 


2 


4 


*di 


2 


1 


1+ 





1/2 





5 


3 ^ 


2 


1 


1+ 





3/2 


2 


6 


3 Po 


1 


1 


o- 


1 


1/2 


1 


7 




1 





1- 





1/2 


1 


8 




1 





1- 





3/2 


1 


9 


3 Pi 


1 


1 


1- 


1 


1/2 


1 


10 


3 Pi 


1 


1 


1- 


1 


3/2 


1 


11 


l d 2 


2 





2+ 


1 


3/2 


2 


12 


l d 2 


2 





2+ 


1 


5/2 


2 


13 


? 'd 2 


2 


1 


2+ 





3/2 


2 


14 


3 d 2 


2 


1 


2+ 





5/2 


2 


15 


3 P2 


1 


1 


2~ 


1 


3/2 


1 


16 


3 P2 


1 


1 


2 - 


1 


5/2 


3 


17 


3 /2 


3 


1 


2" 


1 


3/2 


1 


18 


3 /2 


3 


1 


2~ 


1 


5/2 


3 



TABLE II. Quantum numbers of the three-body channels. 
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^ Meshpoints 


7 < 1+ 


7 < 1 


7 < 2 


6 


-8.3318 


-8.1383 


-9.0551 


12 


-7.3571 


-7.1376 


-7.4330 


24 


-7.3150 


-7.0913 


-7.3688 


36 


-7.3156 


-7.0919 


-7.3688 


40 


-7.3156 


-7.0919 


-7.3688 



TABLE III. Triton binding energies (in MeV) with the Paris (EST) potential (56555557). The 
notation (56...) specifies the employed separable representation as explained in Sec. P4 
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E t (MeV) P(S) P(S') P(P) P(D) 



3 < 1 + 


Paris 


(EST) (11) 


-7.451 


90.63 


1.636 


0.042 


7.692 




Paris 


(EST) (34) 


-7.266 


89.88 


1.652 


0.065 


8.402 




Paris 


(EST) (56) 


-7.316 


89.90 


1.634 


0.064 


8.401 




Paris 


( W-matrix) 


-7.300 


90.22 


1.450 


0.064 


8.265 




Paris 




-7.297 


89.88 


1.625 


0.066 


8.428 




Paris- 




-7.310 


89.88 


1.623 


0.066 


8.428 


i<2+ 


Paris 


(EST) (1111) 


-7.464 


90.62 


1.636 


0.042 


7.704 




Paris 


(EST) (3444) 


-7.375 


89.87 


1.618 


0.066 


8.447 




Paris 


(EST) (5644) 


-7.424 


89.89 


1.601 


0.066 


8.444 




Paris 


(EST) (5655) 


-7.426 


89.89 


1.600 


0.066 


8.446 




Paris 


( W-matrix) 


-7.343 


90.21 


1.436 


0.064 


8.288 




Paris 




-7.408 


89.87 


1.591 


0.068 


8.474 


3 < 1 


Paris 


(EST) (11111) 


-7.464 


90.83 


1.468 


0.044 


7.658 




Paris 


(EST) (34333) 


-7.074 


90.28 


1.492 


0.064 


8.167 




Paris 


(EST) (56444) 


-7.093 


90.30 


1.488 


0.062 


8.154 




Paris 


(EST) (56555) 


-7.092 


90.30 


1.488 


0.062 


8.153 




Paris 


( W-matrix) 


-7.150 


90.51 


1.301 


0.066 


8.123 




Paris 




-7.103 


90.28 


1.468 


0.063 


8.193 


3 <2 


Paris 


(EST) (11111111) 


-7.549 


90.61 


1.459 


0.047 


7.879 




Paris 


(EST) (56555557) 


-7.369 


90.14 


1.420 


0.063 


8.379 




Paris 


( W-matrix) 


-7.088 


90.54 


1.366 


0.062 


8.034 




Paris 




-7.378 


90.11 


1.403 


0.064 


8.418 



TABLE IV. Triton wave function components for the Paris potential. The notation (56...) etc. 



specifies the employed separable (EST) representation as explained in Sec. p4 The result for 



Paris-r is taken from Ref. [26]. 
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E t (MeV) P(S) P(S') P(P) P(D) 



3 < 1 + 


Bonn A (EST) (11) 


-8.350 


92.75 


1.415 


0.028 


5.811 




Bonn A (EST) (44) 


-8.347 


92.35 


1.427 


0.034 


6.188 




Bonn A (EST) (56) 


-8.380 


92.31 


1.432 


0.035 


6.220 




Bonn A (W-matrix) 


-8.371 


92.01 


1.385 


0.041 


6.565 




Bonn A 


-8.378 


92.32 


1.426 


0.035 


6.217 


3 <2+ 


Bonn A (EST) (1111) 


-8.360 


92.74 


1.415 


0.027 


5.820 




Bonn A (EST) (4444) 


-8.411 


92.35 


1.411 


0.034 


6.204 




Bonn A (EST) (5644) 


-8.444 


92.31 


1.415 


0.034 


6.236 




Bonn A (W-matrix) 


-8.399 


92.00 


1.380 


0.039 


6.578 




Bonn A 


-8.443 


92.32 


1.411 


0.035 


6.235 


3 < 1 


Bonn A (EST) (11111) 


-8.298 


92.95 


1.248 


0.030 


5.772 




Bonn A (EST) (44444) 


-8.083 


92.72 


1.252 


0.037 


5.995 




Bonn A (EST) (56444) 


-8.115 


92.68 


1.254 


0.037 


6.027 




Bonn A (W-matrix) 


-8.160 


92.36 


1.209 


0.043 


6.391 




Bonn A 


-8.127 


92.69 


1.248 


0.037 


6.029 


3 <2 


Bonn A (EST) (11111111) 


-8.395 


92.81 


1.264 


0.031 


5.895 




Bonn A (EST) (56444445) 


-8.285 


92.59 


1.236 


0.037 


6.135 




Bonn A 


-8.295 


92.59 


1.231 


0.037 


6.138 



TABLE V. Triton wave function components for the Bonn A potential. (56...) etc. specifies 



the employed separable (EST) representation as explained in Sec. Ill 
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£ t (MeV) P(S) P(S') P(P) P(D) 



3 < 1 + 


Bonn B (EST) (11) 


-8.209 


91.95 


1.361 


0.035 


6.659 




Bonn B (EST) (44) 


-8.137 


91.36 


1.369 


0.047 


7.224 




Bonn B (EST) (56) 


-8.170 


91.34 


1.373 


0.048 


7.236 




Bonn B (W-matrix) 


-8.161 


91.21 


1.286 


0.053 


7.448 




Bonn B 


-8.165 


91.35 


1.368 


0.049 


7.235 


3 <2+ 


Bonn B (EST) (1111) 


-8.219 


91.94 


1.360 


0.034 


6.668 




Bonn B (EST) (4444) 


-8.197 


91.35 


1.354 


0.047 


7.245 




Bonn B (EST) (5644) 


-8.230 


91.34 


1.357 


0.048 


7.256 




Bonn B (W-matrix) 


-8.190 


91.21 


1.286 


0.053 


7.448 




Bonn B 


-8.226 


91.34 


1.354 


0.049 


7.257 


3 < 1 


Bonn B (EST) (11111) 


-8.159 


92.14 


1.200 


0.037 


6.620 




Bonn B (EST) (44444) 


-7.855 


91.75 


1.216 


0.048 


6.985 




Bonn B (EST) (56444) 


-7.884 


91.74 


1.218 


0.048 


6.700 




Bonn B (W-matrix) 


-7.926 


91.56 


1.143 


0.054 


7.247 




Bonn B 


-7.899 


91.74 


1.210 


0.049 


7.001 


3 <2 


Bonn B (EST) (11111111) 


-8.282 


91.96 


1.208 


0.039 


6.791 




Bonn B (EST) (56444445) 


-8.088 


91.61 


1.189 


0.048 


7.149 




Bonn B (W-matrix) 


-7.919 


91.56 


1.143 


0.052 


7.249 




Bonn B 


-8.103 


91.62 


1.184 


0.048 


7.152 



TABLE VI. Triton wave function components for the Bonn B potential. (56...) etc. specifies 



the employed separable (EST) representation as explained in Sec. |III 
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FIGURES 




channel 3 channel 4 



¥(M)[fin1 %,q)[fa] 

0.5 r 1 r 




channel 5 channel 6 




FIG. 1. Graphs of the triton wave function obtained with the Paris (EST) potential. The 
channels are defined in Table 0. 
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channel 7 



channel 8 






FIG. 2. - continued. 
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channel 13 



channel 14 






FIG. 3. - continued. 
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channel 1 



channel 2 






FIG. 4. Same as Fig. 1 but for the Bonn A (EST) potential. 
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channel 7 



channel 8 
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channel 13 



channel 14 





channel 17 




channel 18 



-P(p,q)[fm 3 ] 




FIG. 6. - continued. 
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